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Abstract 

The oscillation frequency of B°-B° mixing (AmJ has been measured using 29.1 fb _1 
of data collected with the Belle detector at KEKB. This measurement is made 
through the distributions of the proper decay time difference of B pairs in events 
tagged as same- and opposite- flavor decays. In each event, one B is fully recon- 
structed in a flavor-specific hadronic decay mode, while the flavor of the other is 
extracted through a likelihood calculated from the 6-flavor information carried in 
its final decay products. We obtain Amd = (0.528 ± 0.017 ± 0.011) ps _1 , where the 
first error is statistical and the second error is systematic. 

PACS: 12.15.Hh, 11.30.Er, 13.25.Hw 



1 Introduction 



In the Standard model, B°-B° oscillation occurs through second order weak 
interactions where the dominant contribution has an internal loop that con- 
tains virtual t-quarks. As a result, the oscillation frequency Arrid is related 
to the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements V t d and V th [1]. 
Since \V t b\ is expected to be very close to 1, a precise measurement of Arrid, 
in principle, provides a method to determine \ Vu\- Although the extraction of 
\V t d\ from Arrid is currently hampered by theoretical uncertainties on hadronic 
matrix elements, precision measurements of Arrid could be important for fu- 
ture determinations of \V t d\- In addition, B°-B° mixing is one of the sources 
of CP violation in B° decays and a good knowledge of Am d is important for 
precise measurements of these CP asymmetries. Since the first observation of 
B°-B° mixing [2], a number of measurements of Arrid have been reported [3]. 

1 on leave from Nova Gorica Polytechnic, Nova Gorica, Slovenia 
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Recently, asymmetric £>-factory experiments operated at the T(4S I ) resonance 
have significantly improved the precision of Am d [4]. 

In this letter, we present a determination of Am d from the time evolution 
of opposite-flavor (OF; B°W) and same-flavor (SF; B°B°, WW) neutral B 
decays at the T(4S f ) resonance. In the absence of background and detector 
effects, the time- dependent probabilities of observing OF (V OF ) and SF (V SF ) 
states are given by 



V OF (At) = J- exp [1 + COS (Am d At)] , (1) 

V SF {At) = exp [1 - cos(Am d At)] , (2) 

4r B o \ t b o J 

where tbo is the B° lifetime, and At is the proper time difference between the 
two B meson decays. As the decay width difference of the two mass eigenstates 
of the B°B° system is expected to be very small in the Standard Model [5], 
we assume in this analysis that it is equal to zero. 

The analysis described here is based on a 29.1 fb _1 data sample, which contains 
31.3 x 10 6 BB pairs, collected with the Belle detector [6] at the asymmetric- 
energy KEKB storage ring [7]. KEKB collides an 8.0 GeV electron beam and 
a 3.5 GeV positron beam with a crossing angle of 22 mrad, resulting in a 
center-of-mass system (cms) moving nearly along the z axis (defined as anti- 
parallel to the positron beam) with a Lorentz boost of (Prfr = 0.425. Since B 
mesons are nearly at rest in the T(4S) cms, the proper time difference At can 
be approximated as At ~ Az/[c(/3^y)r], where Az is the distance between the 
decay vertices of the two B mesons in z (typically 200 fjm). The measurement 
involves the reconstruction of the decay of a neutral B meson in a flavor- 
specific hadronic mode, the determination of the 6-flavor of the accompanying 
(tagging) B meson, the reconstruction of the two decay vertices to determine 
At, and a fit of the At distribution taking into account At resolution and 
backgrounds to extract Am d . We use the same flavor-tagging method as for 
the sin 2<pi measurement [8], and the same vertex reconstruction and resolution 
parameters as those used in the lifetime measurements [9]. 

The Belle detector consists of a three-layer silicon vertex detector (SVD), a 
50-layer central drift chamber (CDC), an array of aerogel Cherenkov coun- 
ters (ACC), time-of-flight scintillation counters (TOF), an electromagnetic 
calorimeter containing CsI(Tl) crystals (ECL), and 14 layers of 4.7 cm thick 
iron plates interleaved with a system of resistive plate counters (KLM). All 
subdetectors except the KLM are located inside a 3.4 m diameter super- 
conducting solenoid which provides a 1.5 T magnetic field. The impact pa- 
rameter resolutions for charged tracks are measured to be a 2 = (19) 2 + 
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(5O/(p/3sin 3 / 2 0)) 2 /im 2 in the plane perpendicular to the z axis and a 2 = 
(36) 2 + (42/ {pf3 sin 5 / 2 0)) 2 /mi 2 along the z axis, where f3 = pc/E, p and E are 
the momentum (GeV/c) and energy (GeV) of the particle, and 6 is the polar 
angle from the z axis. 



2 Event Selection and Reconstruction 

B° mesons are fully reconstructed in the decay modesQ B° — > D + tt~ , D* + tt~ , 
and D* + p~ . Charged pion and kaon candidates are required to satisfy selec- 
tion criteria based on particle-identification likelihood functions derived from 
specific ionization (dE/dx) in the CDC, time of flight, and the response of 
the ACC. Photon candidates are defined as isolated ECL clusters with energy 
more than 30 MeV that are not matched to any charged track. We recon- 
struct 7T° candidates from pairs of photon candidates with invariant masses 
between 124 and 146 MeV/c 2 . A mass-constrained fit is performed to improve 
the 7T° momentum resolution. A minimum tt° momentum of 0.2 GeV/c is re- 
quired. We select p~ candidates as ir~ir pairs having invariant masses within 
±150 MeV/c 2 of the nominal p~ mass. 

Neutral and charged D candidates are reconstructed in the following chan- 
nels: D° — > K~tt + , A"~7r + 7T°, A~7T + 7r + 7r~, and D + — > K~7t + tt + . Candidate 
D* + — > D°7i + decays are formed by combining a D° candidate with a slow 
and positively charged track, for which no particle identification is required. 
We apply mode-dependent requirements on the reconstructed D mass (rang- 
ing from ±15 to ±58 MeV/c 2 ) and the mass difference between D* + and 
D° (ranging from ±3 to ±12 MeV/c 2 ). To reduce continuum background, a 
mode-dependent selection is applied based on the ratio of the second to ze- 
roth Fox- Wolfram moments [10] and the angle between the thrust axes of the 
reconstructed and associated B mesons. 

We identify B decays based on requirements on the energy difference AE = 
£?™- E^ m and the beam-energy constrained mass M^ c = J (E^ m ) 2 — (p ™) 2 , 
where E^ s and p c j^ s are the cms energy and momentum of the fully recon- 
structed B candidate, and E^ m is the cms beam energy. If more than one 
fully reconstructed B candidate is found in the same event, the one with the 
best combined \ 2 f° r AE, M^ c , and the invariant mass of the D candidate is 
chosen. For each channel, a rectangular signal region is defined in the AE-M^ C 
plane, corresponding to ±3cr windows centered on the expected means for AE 
and Mb c . The AE resolution is a = 10 ~ 30 MeV, depending on the decay 
mode and a ~ 3 MeV/c 2 for M bc for all modes (5.27 < M bc < 5.29 GeV/c 2 is 
used). For the determination of background parameters, we use candidates in 

2 The inclusion of the charge conjugate decays is implied throughout this letter. 
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a background dominated (sideband) region, —0.1 (—0.05 for D* + p ) < AE < 
0.2 GeV and 5.20 < M bc < 5.29 GeV/c 2 , and exclude the signal region. 

Charged leptons, pions, and kaons that are not associated with the recon- 
structed hadronic decay are used to identify the flavor of the accompanying 
B meson. We apply the same method that has been used for our measure- 
ment of sin 20i [8]. Initially, the 6-flavor determination is performed at the 
track level. Several categories of well measured tracks that have a charge cor- 
related with the b flavor are selected: high momentum leptons from b — > ci~V, 
lower momentum leptons from c — > s£ + u, charged kaons and A baryons from 
b — > c — > s, high momentum pions originating from decays of the type B° — > 
Jj(*)~X (where X = 7r + ,p + , a*, etc.), and slow pions from D*~ — > D°n~ . We 
use Monte Carlo (MC) simulations to determine a category-dependent vari- 
able that indicates whether a track originates from a B° or B°. The values of 
this variable range from —1 for a reliably identified B° to +1 for a reliably 
identified B° and depend on the tagging particle's charge, cms momentum, 
polar angle and particle-identification probability, as well as other kinematic 
and event shape quantities. The results from the separate track categories 
are then combined to take into account correlations in the case of multiple 
track-level tags. We use two parameters, q and r, to represent the tagging 
information. The first, q, corresponds to the sign of the b quark charge of the 
tag-side B meson where q = +1 for b and hence B°, and q — — 1 for b and B°. 
The parameter r is an event-by-event, MC-determined flavor-tagging dilution 
factor that ranges from r = for no flavor discrimination to r = 1 for unam- 
biguous flavor assignment. It is used only to sort data into six intervals of r 
(boundaries at 0.25, 0.5, 0.625, 0.75, 0.875), according to flavor purity. More 
than 99.5% of the events are assigned a non-zero value of r. 

The decay vertices of the two B mesons in each event are fitted using tracks 
that have at least one SVD hit in the r-<f> plane and at least two SVD hits in the 
r-z plane, under the constraint that they are consistent with the interaction 
point (IP) profile, smeared in the r-(f> plane by 21 /im to account for the 
transverse B decay length. The IP profile is represented by a three-dimensional 
Gaussian, the parameters of which are determined in each run (every 60,000 
events in case of the mean position) using hadronic events. The size of the IP 
region is typically a x ~ 100 /xm, a y ~ 5 jum, and a z ~ 3 mm, where x and y 
denote horizontal and vertical directions, respectively. 

The decay point of the reconstructed B° is obtained from the vertex position 
and momentum vector of the reconstructed D meson and a track (other than 
the slow 7r + candidate from D* + decay) that originates from the fully recon- 
structed B decay. The decay vertex of the tagging B meson is determined 
from tracks not assigned to the fully reconstructed B meson; however, poorly 
reconstructed tracks (with a longitudinal position error in excess of 500 /xm) 
as well as tracks likely to come from K$ decays (because they either form a 
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K$ mass when combined with an another oppositely charged track or miss the 
fully reconstructed B vertex by more than 500 /im in the r-0 plane) are not 
used. 

The quality of a fitted vertex is assessed only in the z direction (because 
of the tight IP constraint in the transverse plane), using the variable £ = 

(l/2n)J2i [(^after — ^before ) /^before] > where n is the number of tracks used in 
the fit, z before and z* fter are the z positions of each track (at the closest ap- 
proach to the origin) before and after the vertex fit, respectively, and £ before 
is the error of z before . A MC study shows that £ does not depend on the 
B decay length. We require £ < 100 to eliminate poorly reconstructed ver- 
tices. About 3% of the fully reconstructed vertices and 1% of the tagging B 
decay vertices are rejected. The proper-time difference between the fully re- 
constructed and the associated B decays, At = t rec — i tag , is calculated as 
At = (z rec — ^ta g )/[c(/?7)x], where z vec and z tag are the z coordinates of the 
fully-reconstructed and associated B decay vertices, respectively. We reject a 
small fraction (~0.2%) of the events by requiring |At| < 70 ps (~45tbo). 

Figure 1 shows the M bc distribution for all the candidates found in the AE 
signal region after flavor tagging and vertex reconstruction. We find 2269 
D + n~, 2490 D* + n~, and 1901 D* + p~ candidates in the signal region with 
average purities of 86%, 81%, and 70%, respectively. The flavor-tagging in- 
tervals I = 1,2, •• -6 contain 2675, 981, 597, 702, 779, and 926 candidates, 
respectively. 



3 Extraction of Am d 



An unbinned maximum likelihood fit is performed to extract Am^, based 
on a likelihood function defined as L(Amd, Wi, W2, ■ ■ ■ ,Wq) = rii-P OF (Ati) x 
]1, P SF (Atj), where the index i (J) runs over all selected OF(SF) events in the 
signal region. The probability density function (PDF)Q, p OF ( SF ) ; i s expressed 
as 



pOF(SF) (At) = _ /ol) pOF(SF) (At . ^ + /ol/ OF(SF) Poi(At)] 

bkg 



+ (1 - /sig)[(l - /oi)P b ° k " (SF) (At) + /oi/ b ° k F g (SF) Poi(At)], (3) 



where / s i g is a signal purity, described below, and W\ is the wrong tag fraction of 
the flavor-tagging interval I containing event i or j. The signal PDF (P^g^ SF ^) 
is the convolution of the true PDF with the resolution function (i? s i g ). The 



Note that the PDF is normalized as J[P OF (At) + P bF (At)]d(At) = 1. 
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background PDF (-Pbkg^ SF ^) is expressed in a similar way: 

+00 

p OF(SF) (At) = J d ( At ')pO F (SF) (A ^ fc(At _ At /) ? (4) 



where k = sig, bkg. A small number of signal and background events have 
large At values. We account for the contribution from these "outliers" by 
adding a Gaussian component P \(At) with a width and contribution fraction 
determined from our B lifetime analysis [9]. The width of the outlier Gaussian 
is cr i = 36 ps, and the fraction, f oh is 0.06% if both vertices are reconstructed 
from two or more charged tracks and 3.1% if at least one vertex is reconstructed 
from a single track constrained with the IP profile. We assume / i to be the 
same for signal and background as well as for the OF and SF sub-samples. 
fk F is the fraction of OF(SF) events for k (= sig, bkg) and given by 
iV° F(SF) /(iV fc OF + iVf F ), where iV° F(SF) is the number of OF(SF) events in the 
signal region for k from the result of the AE-M^ C distribution fit described 
below. 

The signal purity (/ S i g ) is determined on an event-by-event basis as a function 
of AE and M^ c . The two-dimensional distribution of these variables including 
the sideband region is fitted with the sum of a Gaussian function F sig (AE, M bc ) 
to represent the signal, and a background function F^ g (AE, M^ c ), represented 
as an ARGUS background function [11] in M^ c and a first-order polynomial 
in AE. The fraction of signal is taken to be f S i g (AE, Mt, c ) = F S i g /[F s[g + -Fbkg]- 
The fits are done separately for the three B decay modes, and the dependence 
of the signal fraction on the flavor-tag interval I is included in the overall 
normalization of f s \ g . 

In order to account for wrong tagging, the true signal PDF (V^^^) is given 
by replacing cos(ArridAt) with (1 — 2wi) cos(ArridAt) in Eq. 1 (Eq. 2)0. The 
resolution function of the signal is constructed by convolving four different 
contributions: the detector resolutions on z rec and z tag , the smearing of z ta _ g due 
to the inclusion of tracks which do not originate from the associated B vertex, 
mostly due to charm and K$ decays, and the kinematic approximation that the 
B mesons are at rest in the cms. The resolution is determined on an event- 
by-event basis, using the estimated uncertainties on the vertex z positions 
determined from the vertex fit. The parameterization of the resolution depends 
on whether the vertices are reconstructed with multiple tracks or a single 
track. A detailed description of the resolution parameterization can be found 
in Ref. [9]. We use the same parameters for i? s ; g (At) obtained there for this 



We neglect the si^Am^Ai) term that may arise due to the interference between 
Cabibbo- favored and suppressed decay amplitudes [12], since its effect is expected 
to be very small. 
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analysis. The average At resolution for the signal is ~1.56 ps (rms). 



The background PDF is modeled as a sum of exponential and prompt com- 
ponents, 



T> OF(SF)/ A v _ ,OF(SF) 

'bkg y^ 1 ) — /bkg 



f 5 



It bkg 



, (5) 



convolved with i?bkg(At), which is parameterized as a sum of two Gaussians. 
Different parameter values are used for i?bk g (At) depending on whether or not 
both vertices are reconstructed with multiple tracks. The parameters for the 
background PDF are determined using the AE-M\, C sideband region for each 
decay mode. A MC study shows that the fraction of prompt component in 
the signal region is smaller (by ~10-50% depending on the decay mode) than 
that in the sideband region. We correct the estimates of f£ kg for this effect. 

In the final fit, we fix t b » to the world average value [13] and determine Am^ 
and Wi (7=1,6). The fit result is listed in Table 1. Separate fits to the D + tt~ , 
D* + 7i~, and D* + p~ decay modes give consistent Am^ values: 0.536 ± 0.027 
ps -1 , 0.543 ± 0.027 ps -1 , and 0.497 ± 0.032 ps -1 , respectively. Figure 2 shows 
the At distributions for OF and SF events with the fitted curves superimposed; 
Fig. 3 shows the asymmetry between OF and SF events, (N OF — N SF ) / (N OF + 
N$ F ), as a function of |At|. 

The systematic errors are summarized in Table 2. The dominant sources are 
the uncertainties in the resolution functions. The fit is repeated after varying 
the parameters determined from the data (MC) by ±l<r (±2<r). The systematic 
error due to the modeling of the resolution function is estimated by compar- 
ing the results with different parameterizations. The systematic error due to 
the IP constraint is estimated by varying (±10 jum) the smearing used to ac- 
count for the transverse B decay length. The IP profile is determined using 
two different methods and we find no difference between the results. Possible 
systematic effects due to the track quality selection of the associated B decay 
vertices are studied by varying each criterion by 10%. The fit quality criterion 
for reconstructed vertices is varied from £ < 50 to £ < 200. We check the 
systematic uncertainty due to outliers by varying the At range to ±40 ps and 
±100 ps, and find a negligibly small effect. AE-M^ C signal regions are varied 
by ±10 MeV for AE and ±3 MeV/c 2 for M^ c . The parameters determining 
/si g are varied by ±ler to estimate the associated systematic error. We study 
background components with a MC sample that includes both BB and con- 
tinuum events. We find no significant peaking background in the signal region 
above the fitted background curve and conclude that the effect of peaking 
background is negligibly small. The systematic error due to the background 
shape is estimated by varying its parameters by their errors. In the nominal 
fit, we do not include any oscillation component in the background; such a 
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component may arise from the B° originated background. We repeat the fit 
with a background PDF including a mixing term, where background param- 
eters are taken from the sideband data. The dependence on the B° lifetime is 
measured by varying the lifetime by ±ler from the world average value. The 
possible bias in the fitting procedure and the effect of SVD alignment error 
are studied with MC samples; we find no bias. The MC statistical error is 
associated as a systematic error for these sources. 



4 Conclusion 

We have presented a new measurement of Am d using 29.1 ftr 1 of data col- 
lected with the Belle detector at the T(AS) energy. An unbinned maximum 
likelihood fit to the distribution of the proper-time difference of a flavor-tagged 
sample with one of the neutral B mesons fully reconstructed in hadronic de- 
cays yields 

Am d = 0.528 ± 0.017(stat) ± O.Oll(syst) ps" 1 . 

This result has similar precision and is consistent with other recent measure- 
ments at asymmetric 5-factories at the T(4S I ) [4], which are achieving higher 
precision than those at higher energies. Additionally, this measurement con- 
firms the validity of the sin 20! measurement performed on the same data 
sample [8] , since it is based on the same flavor-tagging method, and vertexing 
and fitting procedures. 
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Table 1 



Summary of fit result. 


Fit parameter 


Fit value 


Fit parameter 


Fit value 


Am d 


0.528 ±0.017 ps- 1 






Wl 


0.478 ± 0.017 


W2 


0.313 ± 0.027 


w 3 


0.212 ± 0.030 




0.187 ± 0.027 


w 5 


0.088 ± 0.022 


Wq 


0.016 ± 0.013 



Table 2 



Summary of systematic errors. 


Source 


Error (ps 1 ) 


Resolution parameters 


0.008 


Resolution parameterizations 


0.003 


IP constraint 


0.001 


Track selection 


0.002 


Vertex selection 


0.003 


AE-M hc signal box 


0.001 


Signal fraction 


0.001 


Background shape 


0.002 


Mixing in the background 


0.002 


B° lifetime 


0.002 


Fit bias 


0.005 


Total 


0.011 
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M bc (GeV/c 2 ) 

Fig. 1. Beam-energy constrained mass distribution of the fully reconstructed B° 
candidates. The dashed curve shows the background contributions and the solid 
curve shows the sum of signal and background. 




5 10 15 

At (ps) 



Fig. 2. Distributions of At for (a) SF and (b) OF events with the fitted curves 
superimposed. The dashed, dotted, and solid curves show the background, outliers, 
and the sum of backgrounds and signal, respectively. 
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|At| (ps) 



Fig. 3. Time dependence of the asymmetry between OF and SF events. The curve 
shows the result of the Am^ fit. 
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